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Abstract: We present a simple inflationary scenario that can produce arbitrarily large ^ 

spherical underdense or overdense regions embedded in a standard A cold dark matter ^ ^ 

paradigm, which we refer to as bubbles. We analyze the effect such bubbles would have on 

the Cosmic Microwave Background (CMB). For super-horizon sized bubble in the vicinity 

of the last scattering surface, a signal is imprinted onto CMB via a combination of Sach- dD 

Wolfe and an early integrated Sach- Wolfe (ISW) effects. Smaller, sub-horizon sized bubbles 

at lower redshifts (during matter domination and later) can imprint secondary anisotropies 

on the CMB via Rees-Sciama, late-time ISW and Ostriker-Vishniac effects. Our scenario, 

and arguably most similar inflationary models, produce bubbles which are over /underdense 

in potential: in density such bubbles are characterized by having a distinct wall with the 

interior staying at the cosmic mean density. We show that such models can potentially, 

with only moderate fine tuning, explain the cold spot, a non-Gaussian feature identified in 

the Wilkinson Microwave Anisotropy Probe (WMAP) data by several authors. However, 

more detailed comparisons with current and future CMB data are necessary to confirm (or 

rule out) this scenario. 
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1. Introduction 

Standard inflationary theory predicts that the primordial curvature fluctuations are nor- 
mally distributed ||l], ^, ^, ^ ^. At the level of linear perturbation theory this dictates that 
other observables, such as temperature fluctuations in the cosmic microwave background 
(CMB) should also be normally distributed. This has been confirmed to hold in the real 
data with exquisite precision . 

However, a number of interesting features are present in the data that seems to indicate 
a departure from the minimal inflationary scenario. Wilkinson Map Anisotropy Probe 
(WMAP) Q received a special attention from scientific treasure hunters, since it is the 
only available full sky map of the CMB. In addition to the cold spot (discussed in more 
detail below), various authors have found low multipole alignments 11^, 111], |l^, |l3l, north- 



south asymmetries |14, 15| and various other non-Gaussian features found by blind searches 



T^, 17]. However, the statistical significance of the results is difficult to ascertain, due to 



publication bias and a-posteriori nature of some of the claims presented in the literature 



la 19 



This paper is inspired by the so called cold spot, a nearly-circular region of the CMB 
that is allegedly too cold to be generated by standard Gaussian inflationary perturbations. 
We discuss a generic family of inflationary models/scenarios, as a generalized version of 



multi-stream inflation ||2y, 21 1 , which can produce large spherical regions that have an offset 
in the potential with respect to the surrounding universe. We enumerate possible effects 
that such models would have on the cosmic microwave background. 



These models are thus capable of reproducing the cold spot observed in the CMB 
data, if the bubble is either in the line of sight to the surface of last scattering or actually 
embedded in the surface of the last scattering. We structure the paper as follows: In 
Section Q we discuss the detection of the cold spot and the proposed explanations in the 
literature. In the subsequent section ^ we develop our theory and show that it can fit 
observations without excessive fine tuning (Section ^ and Section |5|). Section ^ concludes 
this article. 



2. Cold spot and proposed explanations 

The cold spot was identified and named in after [^] found a detection of non-Gaussianity 
in the first year WMAP data [p^ j using Spherical Mexican Hat Wavelets [24]. The detection 



significance was about three sigma (about 0.1% probability). The analysis was performed 
on 15 arbitrary wavelet scales and two estimators, skewness and kurtosis, were used. This 
indicates a dilution factor due to a-posteriori detection by a factor of 30 (assuming tests 
are uncor related), lowering the detection threshold to a bit over two sigma. Later analysis 
confirmed this, and estimated that the a-priori detection significance is about 1.9%. pBj . 
The cold spot was also identified using other techniques: using more sophisticated wavelet 
techniques [26|, steerable wavelets [^], needlets ||2^ and scalar indices [^]. On the other 



hand |30| argue that the spot is not statistically significant. Searches were performed look- 
ing for underdensity of astrophysical objects at the location of the spot using radio sources 
and galaxies and found no significant detection [^, |33| . 

The cold spot was found to be almost circular with angular radius of about 4°, tem- 
perature decrement of 70 ixK, and independent of the frequency [p4[ . 

Several different explanations have been proposed for the cold spot. Perhaps the most 
natural is a large void between us and the last scattering surface |35, |3^, which can create 
the cold spot as a secondary anisotropy via the Rees-Sciamma effect. Such void would have 
to be 200-300 Mpc//i in size and have an underdensity of (5 ~ —0.3 at redshift of z ~ 1. 
Such large underdensities are extremely unlikely to appear spontaneously in the standard 
scenario of structure formation. 

Alternative proposals include cosmic texture 35], chaotic post-inflationary pre- 
heating ]39], cosmic bubble collisions (e.g., |4^), and nothing less than a gate to extra 
dimensions ]41]. 
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3. Bubble blowing in mult i- stream inflation 

It is clear that a mechanism that can create nearly circular bubbles or overdensities can 
provide a plausible explanation of the cold spot. To this end we consider a scenario in which 
the inflation can experience different number of e-folds at different spatial positions. This 
can be realized using several different scenarios. Our first example is a two field inflation in 
which the minimum of the inflaton potential causes spontaneous breaking of the mean-fleld 
inflaton path into two distinct paths; an example of such potential is illustrated in Figure 
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Figure 1: Multi-stream inflation. 



|T[ A similar result can be obtained if inflaton is able to quantum tunnel e.g., from stream 
A to stream B in the top panel of Figure ^. 

We can analyze such scenarios without constructing an explicit inflationary model. 
The necessary physics is contained in three parameters: 

• Probability p that the inflaton will bifurcate (or quantum tunnel) into the path B 
rather than remaining on the main path A. For perfectly symmetrical bifurcation in 
the potential p = 0.5, but note that this is not necessarily so. 

• AA^i is the number of e-foldings at which the bifurcation into paths A and B occurs, 
after the start of the observable inflation (i.e. ~ 60 e-foldings prior to the end of 
inflation). 

• is the total difference in the number of e-foldings between the two paths. 

These are parameters illustrated in the Figure 

In what follows, we will denote the regions that follows path B as bubbles, although 
the same mechanism can also produce overdense spherical regions. To estimate the number 
density of such bubbles, we note, that at the time of bifurcation, the number of causally 
connected regions in the volume corresponding to the present-day observable universe is 
given by Hr^e^^^^/H-^, where Hi and Hb are the Hubble constants at the beginning 
of the observable inflation, and the time of bifurcation. Therefore, the total number of 
bubbles following that path B of the inflation potential is given by 

iVobservable = P {Hb/H,fe^^^^ = p e^d-)^^^! = p e'AN* ^ (3_i) 
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Figure 2: Expansion history for multi-stream inflation. 

where we used the definition of the slow-roll parameter e: 

dlnH 



din a 



= {1- e)ANi 



(3.2) 



assuming that e is roughly constant. 

To calculate the physical size of the bubble, we note that the bubble regions reheated 
AN2 e-foldings earlier and did not exponentially expand while reheating. The comoving 
radius of individual bubble regions is thus given by 



-Rbubblo 



1 



-AN^-AN2 



(3.3) 



Finally, we want to calculate the density contrast between the bubble and the rest of 
the universe. In this context, we treat the bubble regions as perturbations in curvature 
and density. Since the curvature perturbation is conserved outside horizon at all orders, 
our calculation holds even for significantly underdense bubbles. 

We start from the separate universe assumption. When we consider super-Hubble 
perturbations, the universe within one Hubble volume is locally an FRW universe, with 
metric 

(3.4) 



ds"^ = AA(x, tfdt^ - a2(t)e-2'/'('''*)dx • dx. 



where Af and ip quantify the local physical time/distance with respect to the temporal 
and spatial coordinates 43|. If perturbations are adiabatic (i.e. pressure is a function 
of density), energy conservations implies that the curvature perturbation of the uniform 
density slices : 

1 dp' 



C(x) =V'(x,t) - - 



(3.5) 



lp,{t) P' + P{p') ' 

is conserved and gauge invariant, where pb{t) and a{t) in Eqs. ( |3.4| )-( |3r5| ) are the background 
(unperturbed) values of density and scale factor. Moreover, p{p) = p/3 after reheating, and 
p = when matter dominates over radiation. One can further show that C equals to the 
e-folding number difference between an uniform density slice and a flat slice. The e-folding 
number in the flat slice is the same as the background e-folding number. Note that the 
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reheating surface is a uniform density slice, so on the reheating surface (the surface with 
temperature T^eh in Fig. we have 



C* = -(iVbubbie - iVo) = A7V2 . (3.6) 

Note that since C* is a constant, the primordial bubble profile is a top-hat. Moreover, is 
gauge invariant. So in a flat slice, we have 

6p = (e-3(i+-)^^2 -l^p, (3.7) 

where w = 1/3 when radiation dominates, and w = when matter dominates. Also, notice 
that when the bubble reheats earlier, we have AiV2 > 0. In this case, 6p < 0. In other 
words, we are producing underdense regions. 

Another important point is that for rare bifurcation events, j3 <C 1, we expect nearly 
spherical bubbles. The reason is that the rare peaks of a random gaussian held with a flxed 
threshold tend to be spherical (e.g., [44|). In our case, the random gaussian field could 



be the transverse fiuctuations about the inflationary trajectory which need to be large to 
trigger bifurcation. 

4. Effects of cosmic bubbles on the CMB 

For the calculation of CMB anisotropics, we shall use the linear metric in the Newtonian 
(or longitudinal) gauge: 

ds'^ = a^{ri) { [1 + 2(/.(x, ri)]drf - [1 - 2V'(x, 7?)](ix • dx} , (4.1) 

where (^(x, rj) and ■0(x, rf) are expressed in terms of comoving coordinates x and conformal 
time ??, which is related to proper time through dt = a{r])dr] (see |]45|| for comprehensive 
review of cosmological perturbation theory). Since (f) = tp + 0{il)'^) for perfect fluids, for 
most of what follows (with the notable exception of Rees-Sciama effect), we use (/> ~ ^. 
The gauge-invariant curvature perturbation: 

is constant on superhorizon scales during the radiation era, and on all scales in the matter 
era, and can be derived from Eq. ( |3.5| ) (for linear perturbations) using the linearized Gqo 
Einstein equation on superhorizon scales: 

-3H{H%b + ijj) = AnGdp, (4.3) 

and Friedmann equation 

H = -AttG{p + p). (4.4) 
Therefore, for the linear Fourier modes V'k relevant for the cold spot we can write: 



Mt) = 9{t)T{\k\)Ck, 



(4.5) 



where T(k) is the transfer function for matter density which goes to 1 on scales bigger than 
the comoving horizon at matter-radiation transition and g{t) is the homogeneous solution 
to the equation of C conservation (Eq. |4.2| ), assuming that g = deep into the radiation 
era: 

2{g + g/H) 
^+ 3(1+^) 

During transition between radiation and matter dominated universes, this equation can be 
solved analytically giving 

/ ^ 3 2 8 2 16 3 / r \ 

where u = r^rad/^^mfl = CLeq/a- Hence, g{t) goes from 2/3 in radiaton era to 3/5 in the 
matter era, and then decays when dark energy starts to dominate. In the transition from 
matter domination to dark energy domination, the equation must be solved numerically. 

4.1 Superhorizon bubbles: Sachs- Wolfe and early ISW effects 

To be detectable in the CMB, the bubble should have a superhorizon size if it is close to 
the last scattering surface. This significantly simplifies the calculations, as the curvature 
perturbation C remains constant. The CMB anisotropics are limited to Sachs- Wolfe and 



early Integrated Sachs- Wolfe (ISW) effects |46]. The latter is due to the fact that close to 



the matter-radiation transition the Newtonian potential is varying. The anisotropics are 
given by: 

^T(A) „ „ ^^oT , gV^(A(r?o-7?),r?) 

__ = 'ilj{rLssn, TiLss) + Qirissn, vlss) + 2 / dr] , (4.8) 

where 0(x, r/) is the intrinsic temperature fiuctuation in the photon field, while r^o and 
r]Lss a-re the conformal times at present and last scattering respectively. Notice that 
partial derivative in Eq. ( [4. 8)) is with respect to conformal time at fixed position, while the 
integral is over the past light cone, which is why it cannot be taken trivially. 
In order to find 0, we can combine the adibaticity conditions: 

5pra = 3/>m0, ^Prad = '^Prad.Q-, (4.9) 

with the superhorizon Gqo equation (^), which yields: 



using Eq. ([4.2| ). Plugging this in Eq. ( [4.8|) , and using: 

d _ d ^dr ^ ^ (4 11) 

dr] dr] dr] ' 

to replace partial derivative with the total (comoving) derivative, we find: 



T 



^LSS 



((^LSSn) + 2 / dr h-\/'tlj{rh,r]Q - r). (4.12) 







Figure 3: Schematic picture of two bubbles tliat intersect tlie last scattering surface. Solid and 
dashed lines show two typical lines of sight that intersect these bubbles. tjlss, Vem and rjex are 
conformal times at last scattering, as well as the entry and exit intersections of the line of sight 
with the bubble. 



Given that on superhorizon scales, V'0 is non-vanishing only at the boundaries of the 
bubble, we can find a closed form for the Sachs- Wolfe -|- early ISW effects: 

-5T(ny 

- ^ J SW+e-ISW 

C* [-0{riLSS - Ven)0{Vex " VlSs) + '^9{Vex)0{r]ex - VLSs) " 'igir]en)0{r]en " VlSs)] , 

(4.13) 

where r/en and r/ex are the conformal times at the moments that photon enters and exits 
the bubble respectively (see Figure while 6 is the Heaviside step function. Deep into 
the matter era, this result reduces to the standard C*/5 for a bubble that intersects the 
last scattering surface, and vanishes otherwise. 

In Figure we plot the profile of such bubble, although note that the edges of this 
profile will be smoothed on the scale of the sound horizon at last scattering, which is a 
about a degree. In order to find this profile, we find rj^n and r]ex for the line of sight at 
angle 9 from the center of the bubble (Figure P), use Eq. ([4.7D to find g{riex) and g{r]en), 
and plug into Eq. (|4l^ ) to find 6T/T. 
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Figure 4: This plot shows the profile of a super-horizon bubble of 1 Gpc/h radius, centered at 
z — 500, 800, 1200 (solid, dashed, dotted) lines. At small distances from the center of the bubble, 
the value of g at the photon exit is constant and corresponds to deep matter era. At the lines of 
sight at the edge of the bubble, g is increasing. At the very edge of the bubble, photons originate 
outside the bubble for low redshift bubbles and hence the profile reflects the difference in the value 
of g between bubble entry and exit. The edges of this profile will be smoothed on the scale of the 
sound horizon at last scattering, which is a about a degree. 

4.2 Subhorizon bubbles: Late-time ISW and Rees-Sciama effects 

For bubbles that had sub-horizon size during the radiation dominated era, the boundaries 
will not be very sharp as the primordial top-hat profile should be convolved with the 
transfer function T(k). Assuming linear evolution, and using Eq. ([4.5|), this implies: 



5(|x|) 



-^bubble 



(2^ 



V^(r,ry) = <7(t)C*5(|r-re|), 
exp(ik-x)r(|k|)Ty(|k|i?bubbie), 



where 



W{x) = ?>x '^(sin X — cos x) 



(4.14) 
(4.15) 

(4.16) 



is the top-hat filter in the Fourier space. We illustrate this in the Figure ^. 



To find the linear late-time ISW effect |46|, it is sufficient to assume the bubble profile 
does not change in a light crossing time. Within this approximation, the impact on the 
CMB just depends on the projection of the bubble profile S{r): 



5T{b) 



^ J usw 
87ra5rC*^bubbic 



dy S (^^2T^) 



d^k 

(2^ 



exp(ik-b)T(|k|)W^(|k|i2b„bbie), 



(4.17) 
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Figure 5: This figure sliows tlie bubble profile in potential (left panel) and overdensity (right 
panel) for C, = 5 x 10"^ for buble size of 100 Mpc/h (solid), 300 Mpc/h (dashed) and 500 Mpc/h 
(dotted) at redshift z = 1. Note that bubbles considered in this paper have the same mean density 
inside the bubble as the mean cosmic density. 



where b is the distance by which the line of sight misses the center of the bubble. 

Since the late-time (linear) ISW effect depends on g, it only becomes significant when 
dark energy becomes important. Therefore, since dark energy is sub-dominant at high 
redshifts, an over/underdensity in the matter era can only contribute to ISW effect at the 
non-linear level. This is known as the Rees-Sciama effect [^ ]. 

In order to find the Rees-Sciama effect, we can use the Gij Einstein equations, deep in 
the matter era, which take the form: 



1 



-AttGT^ =4TTGpma^u^u^, (4.18) 



where ' = dr^, while n*'s are the components of the peculiar velocity, and can be fixed using 
the matter-era Gqh constraint: 



—Hip^i = AnGprnCLU^ ■ 
For spherically symmetric scalar fields X, 



X, 



d fdX dr 
dx^ V dr dx^ 



r dr 



dX 
dr 



and hence Eq. ( [4.19| ) can be used to write the off-diagonal part of ( [4.18 ) as: 



d_ 

dr 



.i9(^-V0 
dr 



while its trace takes the form 



3(^" + 3ai?V') + ^"'|- 
dr 



dr 



(4.19) 



(4.20) 



(4.21) 



(4.22) 



Equation (|]2T|) can be integrated using 1st order solution for ^ and inserted into (|4.22| ), 
yielding an equation for time evolution of ip' at second order. For fixed g{t) = 3/5 in the 



matter era, this can be solved to give 
where 



18C^ 
175aH 



U{\r 



U{r) 



'dS{r)' 


^ 4 


/■~ dr' 


'dSir'Y 


dr 




Ir ~ 


dr' 



(4.23) 
(4.24) 



and we have ignored the transient homogeneous solution to Eq. (|]22l). Similar to ISW 
effect, the Rees-Sciama effect is the line of sight integral of tp': 



mb)] 




T 


RS J 



drj tp' 



1(2) 



36C^ 



r dy U (^/^T^) , (4.25) 

J — oo 



where we have again assumed that the bubble light crossing time is much shorter than the 
Hubble time. As expected, the total effect is dominated by the ISW contribution at low 
redshifts and by the Rees-Sciama at higher redshifts. 

4.3 Subhorizon bubbles: Ostriker-Vishniac effect 

Another secondary anisotropy that could be produced by large bubbles is the Ostriker- 
Vishniac (or the kinetic Sunyaev-Zel'dovich effect) [^]. This effect is caused by the Doppler 
shift of the CMB photons, scattered by free electrons in the universe. It is given by: 



ov 



rLSS ^-j- 

dr — [1 + 5e{rh, r/o — r)] [h ■ u(rn, rjQ — r)] 
* 



(4.26) 



where Se is the over density of free electrons, and dr/dr is the cosmic mean differential 
optical depth for Thomson scattering: 



dr = Ueaaxdr. 



(4.27) 



Here, fie is the mean physical number density of free electrons in the universe, while or is 
Thomson cross-section. 

On the scales that have not gone through shell-crossing, and assuming a uniform ion- 
ization ratio, the electron and total matter overdensities are the same: 6^ — ^ni' Therefore, 
Goi Einstein equation: 

(«5)'., _ ,_^^3, 



9 



(4.28) 



can be used to fix the kernel in Eq. ( [4.261) . After integration by part, this yields: 



ST{h) 



h{l + X)aT 



ov 



SirGrrir, 



I drg-'^ij{rh,ri) — 
Jo drj 



Xe djag) 
drj 



/b(l + X)C.aT [ d 
SnGmp \ drj 



Xe djag) 
a? drj 



] J" dyS [V^T^ 



3/b(l + XX^arHa^ ( d 



IQirGrrir, 



^ [xeH{l -g){l + w)] ] dyS (V^T^) , (4.29) 
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Figure 6: This figure shows the relative amphtudes of Ostriker-Vishiniac and ISW effects as a 
function of redshift. We assume X = 0.75 for hydrogen abundance, and fully ionized medium. We 
ignore the fact that Helium was most likely neutral at redshift bigger than ^ 3. 

where rup is the proton mass, X ~ 0.75 is the Hydrogen mass fraction in the intragalactic 
medium, and Xg is the ionized fraction (assuming both hydrogen and helium are ionized 
equally) and fh = Qh/^m is the baryon fraction. H and w are Hubble constant and total 
equation of state (total pressure divided by total density), respectively. 

Since both ISW and OV effects scale as C* f S, it is easiest to asses the relative strength 
of the two effect by plotting the ratio of relevant quantities. We do this in the Figure |6|. 
We note that at redshifts of 3, the two effects become comparable and that for larger 
redshifts, the OV begins to dominate. 

We plot the individual effects discussed here for three representative cases in the Figure 
1^ The C perturbations were chosen to give a perturbation of 6T/T 10~^ which is the 
same order of magnitude as the cold-spot perturbation. There are several interesting 
aspects to note regarding these plots. First, the ISW effect dominates at small distances, 
but never completely overtakes the RS effect for sensible values of C that can give the right 
magnitude. At the same time, the sign of C determines whether the two effects are going 
to interfere constructively or destructively. Second, the ISW and OV effects dominate at 
small and large redshift respectively and conspire to give similar sized effects. The result 
is that the total effect has a non-trivial angular structure and redshift dependence and this 
work provides clear templates that can be tested against real CMB data. 

It is interesting to note that a sharp change in the ionization fraction Xg, e.g. at 
reionization, could also lead to a sig nificant OV effect in Eq. (p^ ), as was first pointed 
out in lEc 



Note that, for the purpose of calculation of RS effect in Figure |^, we assumed the 





Figure 7: This figure shows the bubble profile associated with ISW (topmost row), Rees-Sciama 
(second row), Ostriker-Vishiniac (third row) and total contribution (bottom row). We show the 
profile for bubbles size 100, 300 and 500 Mpc/h at rcdshift z = 0.5 with = 10""^ (left panels) 
and (t, = — 10~^ (middle panels) and at redshift z = 10 with = 10"'^ (right panels). 



Universe is still in the matter domination at z = 0.5. Given that matter domination 
ends around z ~ 0.5, we expect this assumption to introduce 0(50%) error in RS, and 
0(10 — 20%) error in total temperature profile. 



5. Comparison with Cold Spot 

The mechanism proposed in this section can, in principle, produce bubble regions that 
can be larger or smaller than the horizon. Therefore, we can explain the bubble using 
two different scenarios: bubble between us and the surface of the last scattering and a 
super-horizon overdensity at the actual surface of the last scattering. 
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If a bubble close to or intersecting the surface of last scattering is responsible for 
the cold spot, the 4 degree angular radius of the cold spot requires a comoving radius of 
-Rbubbie 1 Gpc. Figure (Q) shows that a temperature decrement of 70 iJ^K = 3 x 10~^Tcmb 
requires = — 1 x 10"'*. Furthermore, given that there must be one bubble within a Gpc 
of the last scattering surface (at radius of 10 Gpc), one expect A'^observabie ~ 4 bubbles 
within our observable horizon. Plugging these into inflationary scenario of Sec. gives: 



P 



0.15, AiVj* ~ 1, AN2 ~ -1 X 10"^, for superhorizon bubbles. (5.1) 



For bubbles at progressively lower redshifts, the radius of bubble should be smaller to 
match the observed extent of the cold spot. Given the scarcity of anomalies such as cold 
spot, we do not expect the bubble to be closer to us than z ~ 1, assuming that we sit at a 
random point in the universe. This gives a minimum radius of ~ 200 Mpc for the bubble, 
which implies: 

1 < AiVr < 4, (5.2) 



from Eq. (|3J), given that AN2 = C* < 1- 

Interestingly, we can also estimate p, simply based on the angular size of the bubble. 
To do this, we first note that the angular size of the cold spot implies a comoving radius 
of: 

^bubble ~ (1 Gpc)(r/rL5s), (5.3) 

if the bubble sits at the comoving distance r from us. However, if the closest (and thus 
easiest to see) bubble is at distance r, there should be 

A'observable ~ {rLSs/r f, (5.4) 

bubbles within our observable horizon. Now, combining Eqs. ( p.l[ ) and (3^) yields: 

P ~ A^obsorvablc (^^0-Rbubblo)^ ~ 0.04, (5.5) 

given that Hq^ = 3.0 Gpc//i. In other words, bifurcation into path B is roughly a 2a event 
in the inflationary history. 

Figures (0) and (|6|) show that in order to match a temperature decrement of ~ 10~^ 
for bubbles at 0.5 < z < 10: 

C* = AiV2 ~ 10"^ (5.6) 

i.e. the bifurcated inflationary path (path B) expanded 0.1% more that the regular infla- 
tionary trajectory (path A). 

Finally, we should note that in what came above we have ignored the impact of regular 
inflationary gaussian fluctuations on our predictions. The fact is that both Rees-Sciama 
and Ostriker-Vishniac effects include 2nd order cross-terms with product bubble profile and 
standard gaussian fluctuations. While these terms have zero mean, they lead to a random 
component in the bubble profile, which can be potentially seen in higher order statistics of 
the CMB sky (e.g., see |^, |5^). Such treatment is clearly beyond the scope of the current 
paper, but will be possible given the formulae provided in Sec. (Q). 



6. Conclusions and Future Prospects 



In this paper we have introduced a novel mechanism that could lead to over /underdense 
spherical regions/bubbles in the universe, through a probabilistic bifurcation in the infla- 
tionary trajectory. We have then worked out predictions for the imprint of such bubbles on 
the CMB, and estimated approximate values of the parameters of the inflationary scenario 
that could lead to the observed WMAP cold spot. 

If our mechanism is correct and the cold-spot observed in the WMAP data is indeed 
a massive high redshift bubble, there are several observational probes that can test this 
assumption: 

Of course, the most immediate test will be to investigate how well the CMB profiles 
of Sec. @ can fit the WMAP cold spot. Another promising direction is to reconstruct the 
bubble profile by detecting the lensing of the CMB by the large bubble 1 52 1 , and compare 
it to our prediction for density profile (Fig. |5[). Such tests are well within realms of 
experiments such as Atacama Cosmology and South Pole Telescopes |5^, 54 1. Other tests 
would include observing changes in the structure formation within such bubble [55|. In fact, 
using NVSS radio sources, it has been claimed that there is an anomalous underdensity 
of radio sources coincidental with the position of the cold spot in the CMB , but this 
claim has been discarded upon closer inspection However, Figure |5| shows that our 
scenario predicts a compensated bubble wall, rather than a void, which provides a new 
(and unique) template for comparison with galaxy surveys. 
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